Theory of Disordered Itinerant Ferromagnets I: Metallic Phase 
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A comprehensive theory for electronic transport in itinerant ferromagnets is developed. We first 
show that the Q-field theory used previously to describe a disordered Fermi liquid also has a saddle- 
point solution that describes a ferromagnet in a disordered Stoner approximation. We calculate 
transport coefficients and thermodynamic susceptibilities by expanding about the saddle point 
to Gaussian order. At this level, the theory generalizes previous RPA-type theories by includ- 
ing quenched disorder. We then study soft-mode effects in the ferromagnetic state in a one-loop 
approximation. In three-dimensions, we find that the spin waves induce a square-root frequency 
dependence of the conductivity, but not of the density of states, that is qualitatively the same as 
the usual weak-localization effect induced by the diffusive soft modes. In contrast to the weak- 
localization anomaly, this effect persists also at nonzero temperatures. In two-dimensions, however, 
the spin waves do not lead to a logarithmic frequency dependence. This explains experimental 
observations in thin ferromagnetic films, and it provides a basis for the construction of a simple 
effective field theory for the transition from a ferromagnetic metal to a ferromagnetic insulator. 

PACS numbers: 75.20.En; 75.30.-m; 75.30.Ds; 72.15.Rn 



I. INTRODUCTION 

The theoretical treatment of many-fermions systems is 
a hard problem, especially in the presence of quenched 
disorder. Traditional [-.approaches have been Lan- 
dau Fexmi-liquid theorytl and many-body perturbation 
theory.Hu Impressive progress has been made within the 
framework of the latter. The random phase approxi- 
mation (RPA) has been developed for the Fermi liquid 
phase, and similar theories have been used to describe 
magnetic and superconducting phases. In a more re- 
cent development, the corrections to Fermi-liquid the- 
ory known as "weak-localization effects"u have been de- 
rived using many-body diagrammatic techniques Jj Even 
more recently, field-theoretic methods have been applied 
to the many-fermion problem, which have certain advan- 
tages over the traditional techniques. Most importantly, 
they allow for a straightforward application of the renor- 
malization group (RG), implementing, inter alia, an old 
program of describing the various phases o£,many-body 
systems in terms of stable RG fixed points. U So far this 
program has been carried out for clean and disordered 
Fermi liquids.Q It has also been shown that the weak- 
localization effects can be understood as the leading cor- 
rections to scaling near a disordered Fermi-liquid fixed 
point E A further advantage of the field theoretic ap- 
proach is that it allows for a systematic identification 
and analysis of the soft or massless modes that are re- 
sponsible for the long-distance, long-time properties of 
the system. If desirable, it also allows for the derivation 



of effective theories that keep only the soft modes explic- 
itly, while integrating out all other degrees of freedom 
in some simple approximation. This ability to focus on 
soft modes is particularly important for studies of the 
zero-temperature (T = 0) properties of many-body sys- 
tems, since it turns out that there are many soft modes 
at T = that acquire a mass at nonzero temperature. 
Conservation laws, and analogies to classical fluids, are 
therefore of limited value in identifying the soft modes of 
a quantum system at T = 0. 

An example of soft modes that have no analogs in clas- 
sical physics are the diffusive modes that cause the weak- 
localization effects mentioned above, which are long- 
wavelength/low- frequency nonanalyticities in the fre- 
quency and wavenumber dependence of transport coef- 
ficients and thermodynamic quantities. These nonana- 
lyticities, which generically take the form of power laws, 
are examples of a more general phenomenon known as 
generic scale invariance. Correlation functions which, as 
functions of space and time, are power laws and thus ho- 
mogeneous functions (as opposed to, e.g., exponentials), 
contain no intrinsic length or time scales, hence the name 
"scale invariance" . The most widely known example of 
scale invariance, which is caused by soft modes, occurs 
at critical points as a result of the critical modes. Crit- 
ical points are exceptional points in a phase diagram, 
and reaching them requires a fine tuning of parameters. 
Generic scale invariance, on the other hand, does not 
require any fine tuning and is due to soft modes that oc- 
cur over large regions of parameter space, like the dif- 
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fusive "diffuson" modes in disordered electron sytems 
that cause the weak-localization effects. However, generic 
scale invariance is in no way restricted to quantum sys- 
tems; it has been discussed primarily in classical sys- 
tems ranging from classical fluids and liquid crystals to 
sandpilesE 

It is the purpose of the present paper to develop a com- 
prehensive field theoretic method for describing disor- 
dered itinerant ferromagnets, which so far have not been 
studied with these techniques, with particular emphasis 
on the consequences of the soft modes in such systems. 
Specifically, we generalize the theory of Ref. |s| to the 
magnetic case. Simple approximations within this theory 
yield results that correspond to disordered Stoner theory 
for the equation of state, and to disordered RPA-type 
approximations for the transport properties. We then 
study generic scale invariance effects on the background 
of ferromagnetism. An interesting question in this con- 
text is the role of the Goldstone modes or spin waves that 
arise from the spontaneously broken spin rotational sym- 
metry, which represent additional soft modes compared 
to a system without magnetic long-range order. Since 
the soft mode structure of ferromagnets is thus differ- 
ent from that of nonmagnetic systems, one also expects 
generic scale invariance phenomena that differ from the 
usual weak-localization nonanalyticities. Indeed, we find 
that the spin waves contribute to the leading nonanalyt- 
icitic frequency dependence of the conductivity in three- 
dimensions, but not to that of the density of states. Also, 
we find that they do not lead to a logarithmic frequency 
dependence of the conductivity in two-dimensions. This 
explains experimental observations that have found the 
frequency anomaly in thin ferromagnetic films to be the 
same as those in nonmagnetic metals in an external mag- 
netic field. 

This paper is organized as follows. In Sec. |l| we recall 
the general Q-matrix field theory for disordered interact- 
ing electrons that was developed in Ref. ||, and we show 
that this theory allows for a saddle-point solution that 
corresponds to a disordered Stoner theory. In Sec. Ill we 



calculate the spin and density susceptibilities, as well as 
the conductivity, in a Gaussian approximation. This pro- 
duces generalizations of well-known results to the case of 
disordered magnets. In Sec. IV we proceed to perform a 



one-loop calculation of the density of states and the con- 
ductivity, and calculate the contribution of the Goldstone 
modes to the leading nonanalyticities. In Sec. |v| we con- 
clude with a discussion of our results. The ferromagnetic- 
metal-to-ferromagnetic-insulator transition will be dis- 
cussed in a second paperj!3 which we will refer to as 
(II). 



II. MATRIX FIELD THEORY 



A. Q-matrix theory for fermions 

We start with a general field theory for electrons. For 
any fermionic system, the partition function can be writ- 
ten as a functional integral pier fermionic (i.e., Grass- 
mann valued) fields and 



Z= / D[$,1>] exp(S[0,0]) 



(2.1a) 



where S is the action. We consider an action that con- 
sists of a free-fermion part So, a part <Sdi s describing the 
interaction of the electrons with quenched disorder, and 
a part Si n t describing the electron-electron interaction, 



S — So + S ( i 



■s b 



(2.1b) 



Each field ip or carries a Matsubara frequency index 
n and a spin index a =1\ 4= +, — ■ Since we will deal 
with the quenched disorder by means of the replica trick, 
each field also carries a replica index a. It is useful to 
introduce a matrix of bilinear products of the fermion 
fields, 



I -V'itV^t 

-0U02T 
0U02T 

y -0i T 02| 



-01 T 021 
--01T-021 



-V>1 T 021 
"0U021 
0U021 
-01 T 021 



= Q 



12 



•01TV , 2T \ 
01J.-02T 
-■0 U 0) 2T 
01T02T ) 

(2.2) 



where all fields arc understood to be taken at position x, 
and 1 = (rti,ai), etc. The matrix elements of B com- 
mute with one another, and are therefore isomorphic to. 
classical or number- valued fields that we denote by Q$3 
This isomorphism maps the adjoint operation on prod- 
ucts of fermion fields, which is denoted above by an over- 
bar, on the complex conjugation of the classical fields. 
We use the isomorphism to constrain B to the classical 
field Q, and exactly rewrite the partition functionO 

Z = J £[0,0] e s ^ J D[Q]S[Q-B] 

'^[0,0] e s ®M J D[Q]D[A] e 1 ^-^ 

D[Q]D[l]e A ^ Q ^ . (2.3) 



We have introduced an auxiliary bosonic matrix field A 
to enforce the functional delta-constraint in the first line 
of Eq. ( |2.3[ ) , and the last line defines the action A. The 
matrix elements of both Q and A are spin-quaternions, 
i.e. elements of Q x Q with Q the quaternion field. From 
Eq. (2.2) we see that expectation values of the Q ma- 
trix elements yield single-particle Green functions, and 
Q-Q correlation functions describe four-fermion correla- 
tion functions. The physical meaning of A is that its 
expectatio n val ue plays the role of a self energy (see Ref. 
| and Sec. |l^ below). 
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It is convenient to expand the 4x4 matrix in Eq. (2.2 
in a spin-quaternion basis, 



Ql 2 (x) = £ ( T r ® s i) rQl2(x) 



(2.4) 



r,i— 



and analogously for A. Here tq = So = I2 is the 2x2 unit 



matrix, and r ? - = 



—Si = —mt 
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(j = 1,2,3), with 0-1,2,; 



the Pauli matrices. In this basis, i — and i = 1,2,3 
describe the spin singlet and the spin tri plet, respec- 
tively. An explicit calculation using Eq. (2.2) reveals 
that r = 0, 3 corresponds to the particle-hole channel 
(i.e., products V>V0, while r = 1,2 describes the particle- 
particle channel (i.e. , products ipif> or iptp). From the 
structure of Eq. ( |2.2| ) one obtains the .following formal 
symmetry properties of the Q matrices,u 



°o 


= (-r°Q 2 i 


(r = 


0,3) 


5 


(2.5a) 


rQ 12 


= R r+ll r Q 2 i 


(r = 


0,3; i 


= 1,2,3) , 


(2.5b) 


°o 

rVl2 


- °o 


(r = 


1,2) 




(2.5c) 


r Ql2 


— — r^21 


(r = 


1,2; i 


= 1,2,3) , 


(2.5d) 




i n a ia2 
— rV-m-l,- 


"2^1 






(2.5e) 



The star in Eq. ( 2.5ej ) denotes complex conjugation. 

For the purposes of the present paper, we will be par- 
ticularly in tere sted i n th e matrix elements qQ and 3Q. 
From Eqs. (|J) and (Q we find 



°Q 12 (x))£^i 2 JE^w^w) . 



(2.6a) 



Q 12 (x)) a «J 12 i ]T a $ 1<7 (x) Vi CT (x)> , (2.6b) 



where (. . .} denotes an average taken with the full ac- 
tion. From these expressions we obtain various observ- 
ables in terms of expectation values of the Q fields, for 
instance the particle number density n, the frequency or 
energy dependent density of states N, with energy mea- 
sured from the chemical potential /1 (or, at T — 0, from 
the Fermi energy ep), and the magnetization M, 



(2.7a) 

(2.7b) 
(2.7c) 



with /l/,b the Bohr magneton. Here and in what fol- 
lows we denote fermionic Matsubara frequencies by uj n — 
2wT(n + 1/2), and bosonic ones by f2„ = 2irTn. We use 
units such that h — = 1. 



By using the delta constraint in Eq. (2.2) to rewrite 
all terms that are quartic in the fermion field in terms of 
Q, we can achieve an integrand that is bilinear in ijj and 



ip. The Grassmannian integral can then be performed 
exactly, and we obtain for the action A 
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A[Q, A] = „4 dis + Ant + I Tr In ( G^ 1 - / \ 
+ J dx tr (A(x)Q(x) 



Here 



Go 1 = -d T + dl/2m c 



(2.8a) 



(2.8b) 



is the inverse free electron Green operator, with d T and <9 X 
derivatives with respect to imaginary time and position, 
respectively, and m e is the electron mass. Tr denotes a 
trace over all degrees of freedom, including the continu- 
ous position variable, while tr is a trace over all those 
discrete indices that are not explicitly shown. For the 
disorder part of the action one findsa 



AaislQ] = 



1 



itNft 



dx tr (Q(x)Y 



(2.9) 



with Np the density of states at the Fermi l evel in saddle- 
point approximation (see Ref. || and Sec. Ill C below), 
and r the single-particle scattering or relaxation timc.Ej 
The electron-electron interaction Ai n t is conveniently de- 
composed into four pieces that describe the interaction 
in the particle-hole and particle-particle spin-singlet and 



spin-tri 
in Sec 



plot channels.! 



For reasons that will be explained 
HI B below, for our purposes we can neglect the 
particle-particle channel. We th us drop r = 1, 2 from our 
spin-quaternion basis, Eq. (2.4). In particular, we write 
the interaction as the sum of the particle-hole spin-singlet 
and triplet terms, 



.4 (s) - 

•^int — 



A int [Q]=A^ 



a; 



(2.10a) 



A 



(t) _ 



/rfx£ ( -ir £ £ 

r— 0,3 ni, n2,m a 

x [tr ((r r ®ao)Q£n 1+ m(*))] 
x [tr ((r r ® So )^ 2 Q + m ,„ 2 (x))] , 



(2.10b) 



*£(-i) r E EE 

r— 0,3 711,712,771 a i=l 

x [tr ((r r ® Si )Q^, ni+m (x))] 

x [tr {(r r ® 3i )Q a n a 2+m ^))} . (2.10c) 



Here r( s ) > and > are the spin-singlet and 
spin-triplet interaction amplitudes, respectively, is 
responsible for producing magnetism. 
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B. The Stoner saddle point 

We now look for a saddle-point solution of the field 
theory derived in the previous subsection. We are inter- 
ested in an itinerant ferromagnet, i.e. a state where both 
the density of states at the Fermi level and the magneti- 
zation are nonzero. We are further looking for a homo- 
geneous state, so we d rop the real space dependence of 
the fields. Equations (2/7) suggest the following ansatz 
for the saddle-point fields, 



..-ill 



S12 [—S r o6i Q iT, ni + S r3 Si3iA ni ] 



(2.11a) 
(2.11b) 



Substituting this into Eqs. (§]| - fl2.10j) , and using the 
saddle-point condition 6A/SQ = SA/SA = 0, we obtain 
the saddle-point equations 



G n 

F n 

En 
A n 

Here 



i 

2V 

i 

2V 



k 



-2i 
2i 

TtN-pT 



G n - 4:iT {s) Tj2> 



i r 

Lx. n 



Fn-Air^T^e 1 "™ F n 



^«(k)=2 K(k) + a ? 7(k)] 
are Green functions given in terms of 

s« ( k ) = ■ r~TA — v" 



(2.12a) 
(2.12b) 
(2.12c) 
(2.12d) 

(2.13a) 
(2.13b) 

(2.13c) 



with = k 2 /2m — [i. For later reference we also define 
a matrix saddle-point Green function 



G sp — {Gq 1 — zA 



(2.13d) 



whose matrix elements are given by 



(G sp )i2 = 6 U [g ni {t 8> s ) + T nx (r 3 <8 s 3 )] . (2.13e) 

Clearly, E and A are self energies that describe the 
interaction in Hartree-Fock approximation, and the dis- 
order in self-consistent Born approximation. The Green 
functions Q and T obey the equations 

(*w„ - 6c - £„) S»(k) + A„ jF„(k) = 1 , (2.14a) 
(»Wn-&-Sn) ^„(k)+A„a„(k)=0 . (2.14b) 



If we ab sorb t he second, Hartree-Fock, contribution to E 
in Eq. ( [2.12c ), which is purely real and frequency inde- 
pendent, into a re definition of the chemical potential, we 
can rewrite Eqs. ( 2.12c , 2. 12d ) in the form 



A n = A - 



1 1 



ttN f t V 4^ 

k 



where 



(2.12c') 
(2.12d') 

(2.15) 



From Eqs. (fjcj), ( |2.1l| ), ( |2.12j ), and Q2.15| ) we see that A 
is related to the magnetization in saddle-point approxi- 
mation, A = M/hb- 

To discuss our saddle-point solution, let us first con- 
sider the c l ean l imit, r = oo, £„ = 0, A„ = A. From 
Eqs. (|2l^ , |2~11 ) we then obtain an expression for [the 
equation of state that is familiar from Stoner theoryjLa 



n k 



CO 2 



(2.16) 



In particular, the threshold value of rw for the onset of 
ferromagnetism is given by the Stoner criterion 



n f r (t) = l 



(2.17) 



and the usual physical interpretation of Stoner the- 
ory applies.t-3 For instance, the (k) are Hartree-Fock 
Green functions with the chemical potential shifted by 
±A n , which is essentially the magnetization. 

The discussion of the disordered case proceeds in exact 
analogy to the BCS-Gorkov theory of disordered super- 
conductors (this is the main reason why we have cho- 
sen to write Stoner theory in a form analogous to BCS 
theory) B The equation of state then takes the form 



1 = _ 2T (t) T yl_ 



fE 



1 



£k+27Sgnw„) -A 2 



(2.18) 



The integral is independent of the disorder to lowest orde r 
inl/r, and so is therefore the Stoner criterion, Eq. (2.17). 
This is the magnetic analog of Anderson's theorem in 
superconductivity. 

We conclude that our saddle-point solution of the field 
theory, Eqs. fl2.g|) - ( 2.10 ), describes a ferromagnetic state 
in a disordered Stoner approximation that is analogous 
to BCS-Gorkov theory for disordered superconductors. 



III. GAUSSIAN APPROXIMATION 
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A. Gaussian action 



We now write Q = Q sp + SQ, and A = A sp + SA, and 



expand the action, Eq. ( 2.8a ), in powers of SQ and <5A. 
To Gaussian order we obtain A — A sp + .Ag with 



Aq = Aiis [SQ] + Ant [SQ] + - Tr ( G sp SA G sp 5 A 
+Tr (SA SQ) . 



(3.1) 



Defining Fourier transforms of the fields, A(k) = 
J c?x exp(ikx) A(x), and analogously for Q, the piece of 
Ag that is quadratic in A can be written 



i7EEEE^A) 12 (k) ;iv 34 (k) i(^) 34 (-k) , 

(3.2a) 



V 



k 1234 rs ij 



with the matrix A given by 



%A 12M (k) = 5 13 S M [ <„ 2 (k) m°° Stij + < na (k) mf Siij 
+vl°n 2 (k) m™ y + ^ n2 (k) m^ ;ij ] 
= 5 13 <52 4 ; j > 12 (k) . (3.2b) 



Here we have defined convolutions of Green functions 

^ k ) = -^E G «(p) G ™(p+ k ) > ( 3 - 3a ) 



where i, j = 0, 3 and G° = G, G 3 = J-. The matrices 



,00 



etc. are defined as 
no 



77? 



m 



m 



03 

rs,ij 
30 

rs, 27 



-tr (rrTar^tr (s 3 s 4 s]) 
itr (r r T3r])tr (s;s 3 4) 



,33 



where | | 



4tr(r r r+) 



tr(sist) 



(3.3b) 
(3.3c) 
(3.3d) 

(3.3e) 



= Sio + 5n + 5i2 — <5.j3 , etc. Q and A can 



now be decoupled by shifting and scaling the A field. If 
we define a new field A by 



<5A(k) = 2A- 1 (SA{k) - SQ(k)) , 



(3.4) 



then A and Q decouple. We can thus integrate out SA 
and obtain the Gaussian action entirely in terms of Q, 

MQ] = ~y E E E ErW)"W iha-'umv 

k 1234 rs ij 

xl(5Q) M (-k) + A dis [5Q] + A int [SQ] , (3.5) 



with A 1 the inverse of the matrix A defined in Eq. 
(|3.2b|) . It is convenient to rewrite this result as 

Ag[Q] - ^ E E E E r W)l2(k) ^-Ml2,34(k) 
k 1234 rs ij 



xi(SQU(-k) , 



(3.6a) 



where 



A^i2,3 4 (k) = S 13 S 2i [(An)- 1 ^) - 1/ttNft] 

^1 — 2,3—4 Voc\ai.2 V(xiQL3 

B , (3.6b) 

with 1 the unit 8x8 matrix, and B a matrix whose 
elements are 

B rStij = -S r .8 ij 2TT® , (3.6c) 

where r<°) = -T^ and rt 1 ' 2 ' 3 ) = 

B. Gaussian propagators 

The Gaussian Q propagators are give n in terms of the 
inverse of the matrix M defined in Eq. ( 3.6b| ). As in Ref. 
|§| we solve the redundancy problem inherent in the matrix 
field theory due to the symmetry relations, Eqs. (pH|), 
by formulating the theory entirely in terms of matrix 
elements Q12 with n\ > n-i- Then we have 



( , r (5Q)i2(k)i(5Q)34(-p)> G =«yk,p^;iii 2 ; j > 12 ' 3 4(k) ) 

(3.7a) 



where 



with 



hi = 1+^12 [-i+4i] 



(3.7b) 



J rs = 7 tr TVTsTj S r0 



tr s^Sj s ^ 



t 



tr SjSjS 



1 



'r's '3 ^r3 

-2 <5 rs <5jj [(5 r0 <5io + <5r3 (1 - &o)] 



tr s jSi s 3 



(3.7c) 



Here and in the following, (. . .)q denotes an average 
taken with the Gaussian action. To determine the inverse 
of M. we no tice that the only nonzero matrix elements 
of l J s A, Eq. ( 3.2b ), are the diagonal and the antidiagonal 
elements in our r ® s basis, and that the same is true for 
the inverse of A. We find 

rs-^12,34(k) = <5l3 S 2 4 reruns ( k ) + <5i_ 2 ,3-4 S aia2 S aiCt3 

(k) , (3.8a) 



x v E 

rs^ nin2^n^ri4 y 



where 
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E, 



111112,113114 



(k) = C nin2 (k) B 



-1 ^ ^ni — T12 ,'^5 — ^6 ^n^ne (k) -B 



Cn 3 n 4 (k) 



(3.8b) 



Here i? is the matrix given by Eq. (3.6c), and C is defined 
as 



C nm (k) = (™ 00 + r» 
+ (m 00 
+ (m 03 + ? 
+ (m 03 - ? 

Here we have defined 



)(D+ m (k) + p- m (k))/4 
m 33 )(£+ m (k)+£ nm (k))/4 
m 30 )(P+ m (k)-I?- m (k))/4 
m 30 )(f+ m (k)-C m (k))/4 



(3.8c) 



V 



± 

n m 

c± 

-'nm 



with 



(k) 
(k) 



^ m (k)/[l-^ m (k)/7r7V F r] 
7?± m (k)/[l-7 7 ± m (k)/ 7 r7V F r] 



in terms of 



± 

nm 



(k) 
(k) 



(k) 
(k) 



V>+ m (k)±V+ m (k) 

Vnm(k) ±^nm( k ) 



^° m (k)±^ 33 m (k) 
^ 3 m (k)±^ 30 m (k) 



(3.8d 
(3.8e 



(3.8f 

(3.8, 



(3.8h 
(3.8i 



Let us discuss these results. First of all, we notice 
that setting the magnetization equal to zero results in 
V + = V~ = £ + = £~ = T>, and we recover the results 
of Ref. |^. In particular, T> nm is diffusive for nm < 0, 
and hence the spin-singlet and spin-triplet channels of 
the Gaussian propagator C are all soft. In the magnetic 
case, this changes. T> + and T>~ are still soft, and given 
by the nonmagnetic result with the Fermi energy shifted 
by ±A, respectively, 



^ m (k) = 



D±k 2 + \n n . 



This holds for nm < 0, and in the limit of small frequen- 
cies and wavenumbers. For nm > 0, T>^ is finite in that 
limit. and are the density of states at the Fermi 
level, and the Boltzmann diffusivity, respectively, of an 
electron system whose Fermi energy has been shifted by 
±A. The spin-singlet, and the longitudinal component 
of the spin-triplet, are thus still diffusive, as one would 
expect. However, the transverse component of the spin- 
triplet is massive. A calculation yields, for nm < 0, 



£nm(k) 



(-^-I)+(|fi n _ m |+£>±k2)2± 



where we have defined 



1 ± 2iAr 



(3.10a) 



(3.10b) 



The mass of these "spin-diffusons" is thus proportional 
to the magnetization. However, the "interacting" part 
E of the propagator, Eqs. (3.8), is still soft even in the 
transverse spin-triplet ch anne ls. To see this, consider the 
saddle-point equations ( |2.14 ). Multiplying the first of 
these equations by J- n {\s) and the second one by C/„(k), 
subtracting the second equation from the first one and 
integrating over the wave vec tor, w e obtain <^ n (k = 0) = 
2iF n /A n with F n from Eq. ( ^.12b| ). Using this, it is easy 
to show that 



T^£,±(k = 0) = -l/2r 



(t) 



(3.11) 



Therefore, l J s E is still massless for i,j = 1,2. These are 
of course the magnetic Goldstone modes, and the above 
calculation just proves that our Gaussian approximation 
is conserving in the sense that it correctly reflects the 
symmetries of the problem, and the resulting soft modes. 
We note that there is no frequency restriction on the 
softness of the Goldstone modes, -E12.34 is m assl es s bot h 
for n xn 2 < and for n x n 2 > 0. From Eqs. (|3~8l) , ( |3.10| ) 



and (3.11) we see that the structure of the Goldstone 
modes is 



9 (P,^™) 



l + 2r(')T£„ 3 ,„ 4 <S„3- n4 ,n^ 3 S?(p) 
d 



±ifL 



(3.12a) 



with c and d constants. The second equality in Eq. 
( |3.12a ) holds in the limit of small wavenumbers and fre- 
quencies. The magnetization and disorder dependence of 
c and d is quite complicated, and we write down only the 
clean limit results to leading order for small values of A, 
where one finds 



A/6k 



d = -A/2iV F r ( 







(3.12b) 



For later reference we also determine the Gaussian A 



(3.9) propagator. Using Eq. ( p.4| ) in ( |3.2a| ), we find 



(*(M) 12 (kK(£A)34(-p)> G = 6 k , p ^%h2 %A 12M 



(k) 



(3.13) 



Notice that t he <p ' J 9 (k), Eq. ( 3.3a| ), and therefore the ma- 
trix A, Eq. (3.2a), reduce to finite numbers in the limit 
of low frequencies and small wavenumbers. The A prop- 
agator is thus massive. 

Finally, we mention that if one keeps the particle- 
particle or Cooper channel, the corresponding propaga- 
tors have a mass proportional to the magnetization, just 
like the Cooperons for nonmagnetic electrons in an exter- 
nal magnectic field are massive.Ej Since we are interested 
in universal phenomena that are due to the soft modes in 
the system, this justifies our having neglected the Cooper 
channel. 



G 



C. Physical correlation functions 



X± (q,O n )=T^©± +niro (q) 



(3.17b) 



We now use the results of the preceding subsections 
to calculate some correlation functions of physical in- 
terest. Let us start with the single-particle density of 
s tates (DOS) in saddle -point approximation. From Eqs. 
( |2.7rj|j2.1l4p.l2a|j2.13c| ) we find for the DOS as a function 
of imaginary frequency 



N(iu n ) = - [N HF (n + A„,iw„) + 7Vhf(m - A n ,iui n )] . 

(3.14) 

Here -/Vhf(amw„) is the DOS for nonmagnetic electrons 
with chemical potential \i in the disordered Hartree-Fock 
approximation of Ref. |8|. In particular, N-p = 7Vhf(Mi 
The DOS as a function of the real frequency is obtained 
as N (u>) = lmN(iuj n — * u> + iO). We recognize in Eq. 
( |3. 14 ) the usual result of Stoner theory. The physical in- 
terpretation is that the band gets split into a band of up- 
spin-electrons and a band of down-spin-electrons, with, 
the energy splitting proportional to the magnetization.Ej 
Also of interest are the density susceptibility, and the 
longitudinal and transverse spin susceptibilities. Let us 
define a general susceptibility 

£x(q,0„)=16T J2 (i(tQ)Z-n,mM) 



xi(SQ) 



a a 

m-2 ,m 2 +n 



(-q)> 



(3.15) 



By adding an appropriate source term, or by writing the 
density susceptibility Xn as a fermionic expectation value 
and translating to a Q-field correlation by means of Eq. 
Q, we find 

Xn (q, n n ) = °° X (q, fin) + §§x(q, On) - < osXfa, fin) 



^oX(q,fin) 



(3.16a) 



Similarly, one obtains the longitudinal (L) and transverse 
(T) spin susceptibilities as 

XL(q, fin) = TOX(q, fin) + 33X(q> fin) ~ *03X(<l>fin) 



«30X(q, fin) 



(3.16b) 



XT (q, fin) = ooX(q, fin) + 33%(q> fi n) ~ % OsXfa, fin) 



«30X(q, fin) 



(3.16c) 



We can calculate these susceptibilities explici tly in 
Gaussian approximation by using the results of Sec. Ill B. 
For the density susceptibility we find 



Xn 



X++X-+ 4TW x+ 



X- 



i + (rw - rw)( x+ + x-) - 4rwr(*) x+ x- 



(3.17a) 



Here all susceptibilities are understood to be functions of 
q and fi„, and the 



with T> from Eq. ( 3.8d ), are disordered Lindhard func- 
tions for nonmagnetic electrons with the chemical poten- 
tial shifted by ±A„. For small frequencies and wavenum- 
bers, and to leading order in the disorder, they read 



X±(q, 



N F D ± q 2 



D±q 2 



with 



D± = ^( M± A) 



(3.17c) 



(3.17d) 



diffusion constants for electrons in the spin-up and spin- 
down bands, respectively, in d spatial dimensions. Note 
that our Gaussian approximation respects particle num- 
ber conservation, as expressed by the fact that Xn(q - ► 
0, Q n ) — > 0. However, the structure of Xn is not diffusive. 
Rather, in the limit of small q and Q n it takes the form 



— jfl„q 2 ai 



q 4 a2 



(-ifi n ) 2 - «S!„q 2 a 3 + q 4 a 4 



(3.18) 



with ai , a2, £13 , and four independent parameters that 
depend on D+ , D~, N F T (s \ and T (t \ 

From Xn we also obtain the conductivity through the 
identity 



a(ifl n ) = lim -^Xn(q, «^„) ■ 
q^o q z 

In the static limit we obtain <r(0) = ctb, where 



D~ 



N V D 



(3.19) 



(3.20) 



with D the Boltzmann diffusion constant of a nonmag- 
netic system with unshifted chemical potential. In Gaus- 
sian approximation, the conductivity of our itinerant fer- 
romagnet has thus the ordinary Boltzmann value. 

The longitudinal spin susceptibility can also be 
expressed in term s of x+ an( i X-- From Eqs. 
(|3T6b| , |3.15| , ^7l jg) we obtain 



XL 



x+ + x 



+ y- -4r( s ) 



x+x- 



i + (rw - rW)(x+ + x-) - 4r«r«x+x- 



(3.21) 



We note again that xl, like Xn, is massless, but not dif- 
fusive. In the clean limit, t — > oo, x+ an d X- become 
Lindhard functions proper with shifted chemical poten- 
tials. pLn that limit we recover the RPA result of Izuyama 
et al. ,E3 if we take into account that their Hubbard model 
with coupling constant v corresponds to the special case 
FW = r^ s ) = — v/2 in our notation. Another special case 
is the nonmagnetic one, where x+ = X- — Xo ; with xo 



7 



the (disordered) Lindhard func tion with chemical poten- 
tial \x. In that case, Eq. ( 3.21 ) reduces to the ordinary 
RPA result, X l = 2 Xo /(l + 2T~Wx ). 

Finally, we calculate the transverse spin susceptibility. 
We find 



x T (q,n„) = 



m m+n,m 



(q) 



i + 2r«r£ ro £+ +n>m (q) 

^Em £m+n,m (q) 



i + 2r(*)TE m £ ; 



m m+n,m 



(q) 



(3.22) 



with £± from Eq. (3.8d) 



To discuss this resul t, we use again Eq. ( |3.1l[) . Insert- 
ing this into Eq. ( |3.22| ) yields l/x T (q = 0, Q n = 0) = 0, 
which is indicative of the Goldstone modes. An expan- 
sion in powers of q and f2 n yields the familiar quadratic 
dispersion relation for magnons in itinerant ferromagnets, 
viz. 



-1 



(3.23) 



with ^ from Eq. (3.12a). In the limit |k|/A:F, ^n/fF << 
A/ep << 1, this resuitagrees with the one obtained with 
elementary methods.Eil 



IV. EFFECTS OF THE GOLDSTONE MODES IN 
PERTURBATION THEORY 

In this section we examine the contributions of the 
ferromagnetic Goldstone modes (FMGM) or spin waves 
discussed in the previous subsection to the DOS, N(u>), 
and the electrical conductivity, cr(f2). In particular, we 
focus on the contributions of these soft modes to the non- 
analytic frequency dependences of these quantities. Our 
goal is to determine whether or not these additional, com- 
pared to a Fermi liquid state, soft modes that are due 
to the long-range ferromagnetic order, contribute terms 
to N(uj) and a (ft) .that are as strong as the usual weak 
localization effects.El We will show they do in the case 
of the conductivity in three-dimensions, but not in two- 
dimensions, and not in the case of the DOS in any di- 
mension. This strongly suggests that the FMGM are not 
important, at least near two-dimensions, in determining 
the properties of the metal insulator transition from a 
ferromagnetic metal to a ferromagnetic insulator. This 
last point will be important in (II). 



A. Single-particle density of states 

The c ompu tation of the DOS is straightforward us- 
ing Eq. ( 2.7b[) a s a starting point. We first write the 



action, Eqs. ( |2.8[ ), as th e sa ddle-point contribution, plus 
the Gaussian part, Eq. (|3.l|), which is quadratic in Q and 



A, plus non-Gaussian (cubic and higher order) terms, 



o 



(a) 



CD 



(b) 



FIG. 1. One- loop (a) and two- loop (b) contributions to the 
one-point vertex function. 



A[Q, A] = A[Q sp , A sp ] + Ag [SQ, SA] +J^A e [S A] 



e=3 



where 



A i [SA] = --Tr(iG sp SA) i 



(4.1a) 



(4.1b) 



By introducing SA again as defined in Eq. (3.4), we can 



decouple SQ and SA in the Gaussian term, at the expense 
of having the higher order terms depend on both Q and 
A, 



A[Q, A] = Asp + Ag[SQ, SA] + ^ Ai[2A^ 1 (SA - SQ)] . 

(4.1c) 



e=3 



We proceed by writing the one-point Q-correlation 
function on the right hand side of Eq. J2.7b ) as, 



(4.2) 



The diagrammatic loop expansion for the irreducible part 
of (SQ), or equivalently for the one-point vertex function, 
is shown in Fig. [j]. The one-loop term, Fig. |l|(a) , is given 
analytically by 

(SQZ(x))i-ioo P = (SQZ(x)A 3 [2A-i(6A - SQ)] > G . 

(4.3) 

Equation ( |4.3| ) is evaluated by usi ng W ick's t heor em and 
the Gaussian propagators, Eqs. (3/7) and ( 3.13 ). The 
leading nonanalyticities are given by contributions where 
the loop in Fig. |l|(a) is a soft mode. That is, the loop 
must be a Q- propa gator, since the A-propagator is mas- 
sive, see Eq. (3.13). As we have seen in Sec. [II, there 
are two different types of soft Q-propagators, viz. the 
diffusive modes or diffusons in the spin-singlet and lon- 
gitudinal spin-triplet channels, and the Goldstone modes 
or spin wave s in the transverse spin-triplet channels. The 
former, Eq. ( |3.9| ), contribute a term that is very similar 
to the weak localization contribution to the DOS for non- 
magnetic electrons in an external magnetic field. We will 
discuss this contribution in detail in (II). For now we 
just consider the functional form of these contributions. 
They take the form of a frequency-momentum integral 
over a diffuson propagator. The most divergent part of 
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the latt er is the interacting piece, which is denoted by E 
in Eq. ( 3.8b ). Inspection shows that this piece is essen- 
tially a diffusion propagator squared, i.e. it scales like 
1/lu 2 ~ 1/k 4 (here, and in similar arguments below, k 
and u> denote generic internal wavevectors and frequen- 
cies.) Dimensional analysis then yields 



N(fl + w)l-loop,diffusons = 0(lO 



(d-2)/2 ) 



(4.4a) 



The other soft mode contribution is due to the FMGM, 
or spin waves. They take the form of a frequency- 
momentum integral over the FMGM propagator, Eq. 
(3.12a), which scales like 1/k 2 ~ We thus have 



N(fi + w)i_loop,FMGM = 0(uJ d/2 ) 



(4.4b) 



We conclude that the FMGM contributions to the DOS 
are subleading compared to the diffuson contributions. 
We will discuss this result further in Sec. ^ below. Here 
we just mention that it implies that the FMGMs are ir- 
relevant, as far as the DOS is concerned, for all leading 
universal effects in the limit of long wavelengths and low 
frequencies, in particular for the critical behavior at the 
metal-insulator transition from a FM metal to a FM in- 
sulator near d = 2. 



B. Electrical conductivity 

We now investigate the same point for the electrical 
conductivity. In terms of Grassmann variables, the dy- 
namical electrical conductivity as a functipn.-af the real 
frequency is given by the Kubo formulaJ13ll3 



cr(O) 



in/m e 

n 



iT 



nvm 2 



x(c iff (x)a xi < +mifl (x)c 2a ,(x') 

^<< 2 -^(x')>L m ^ + io ■ ( 4 - 5 ) 

Here x = (xi, xi, X3), etc. The gradient operators in Eq. 
(4.5) imply that cr(fi) cannot be written as a simple Q- 
correlation function. At this point we have two choices. 
We can either generalize the fields we work with, or, we 
can introduce a suitable source field to generate a. We 
choose the second path, following Refs. pJLgOl 



1. Source formalism 



To the fermionic action, Eq. ( 2.11: ), we add a source 
term Sj, 



s i = - E E^> / dx E E d ^ c^w 



(4.6a) 



After integrating out the Grassmann fields, the action 



A[Q,A,j] is given by Eq. (2.8a) with the Tr In term re- 
placed by, 



Ai n [A,j] = -Trln(Go 1 +JL-iA) , (4.6b) 



with 



where 



(JL) 12 = ^j^L 12<m d xl , (4.6c) 



L\2,m — 0~a 1 a 2 [(t_® So) ^ni,n 2 +m ~~ { T +® s o) ^na.ni+m] i 

(4.6d) 

with 

r± = t ± ir 3 . (4.6e) 



In the presence of the source, the partition function, Eq. 
(2.3), turns into a j-dependent generating functional, 



Z\j] = / D[Q] / D[A] e^«' A J] 



(4.7a) 



In terms of derivatives of Z\j], a(Q) is given by, 
in/m e iT d 2 



a(Sl) = — 



O 



nVm 2 e dj«dj% 



lnZ[j] 



3=0 



(4.7b) 



There are sever al pos sible ways to proceed at this 
point. References [19 20 derived a nonlinear a model. 
Here, we will use instead a si mple p erturbation expan- 
sion. We first expand A\ a , Eq. (4.6b), in powers of j and 
SA: 



An [A] = A (V 



(4.8) 



i,t=0 



with j4(i^) = 0(j l , SA e ). According to Eq. (4.7b) we can 
restrict ourselves to terms with i < 2. The saddle point 
contribution is, 



.4 



(2,0) = ~ Tr (G sp J LG sp J L) 



(4.9) 



In th is ap proximation, and using E qs. (2.13) and ( |4.6| ) in 
Eq. (4.9), and the result in Eqs. (4.7), one obtai ns th e 
Boltzmann equation result for a as given by Eq. ( 3.20| ). 

For the fluctuation corrections to the Boltzmann con- 
ductivity, there are two classes of terms: An,t) an d A^,e) , 
both with £ > 2. Since there are two j-derivatives in 
Eq. ( 4.7b ) for er, the relevant correlation functions are 



(A 2 X £ j)g and (A(2,£))g- The diagrams that give the one- 
loop contributions to a are shown in Fig. |^. Analytically 
one needs, 
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a = o B + 




+ 




+ 0(2-loop) 



FIG. 2. Diagrammatic contributions to the conductivity. 
Wavy lines denote sources JL, and solid lines denote 5 A prop- 
agators. 



-4(1,2) = --Tr(G sp JLG sp 8AG sp 8A) , (4.10a) 



and 



.4 



(2,2) 



Tr (G sp JLG sp 8AG sp JLG sp 8A) 



(4.10b) 



+Tr {G sp JLG sp JLG sp 8AG sp SA) 
Now we write the loop expansion of the conductivity as 

oo 

<j(n^0)=a B +Y J M^^0) , (4.11) 

i=l 

with <7i the i th term in the loop expansion. To one-loop 
order, Eqs. ( f4.7b| ) and ( O ) give, 



<ri(fi) 



a 2 



[o-(l.l) +0"(l,2)] i: 



(^2,2>G + ^ (^(1,2))g 



(4-12) 



where (. . .) c denotes that only connected diagrams are 
to be taken into account. 17(1,1) an d <T(i 2 ) arc defined, in 
order, as the two terms in the previous equality. They are 
graphically represented by the first and second diagram, 
respectively, in Fig. 0. 



2. Analysis of the one-loop terms 

Th e evaluation of the one-loop terms given by Eq. 
(4.12) is straightforward but tedious, crn,i) and 07 li2 ) 

are given by two- and four-5A correlation functions, re- 
spectively. Due to the coupling between 8 A and 8Q these 
correlation functions contain both diffuson and FMGM 
contributions, and hence are soft. 



We first consider <7(i,i). Structurally, this term is anal- 
ogous to the one-l oop c orrection to the DOS that was 
calculated in Sec. [VA. Therefore, while the diffuson 



contribution to this diagram leads to the usual weak- 
localization correction of order f2( d_2 )/ 2 in 2 < d < 4, 
and In 17 in d — 2, the FMGM contribution is weaker, of 
order (Q d / 2 ) in 2 < d < 4. 

The 17(1,2) contribution, which is represented by the 
second diagram in Fig. |^, is more complicated. From 
the spin structure of this diagram it follows that the two 
propagators that form the loop must either both be diffu- 
sons, or both FMGMs; the contribution that mixes these 



two modes vanishes. The diagram with two diffusons 
yields again a weak-localization term of 0(Q( d ~ 2 ^ 2 ), and 
will be discussed in (II). The one with two FMGMs has 
the following structure. Each mode contributes a factor 
of 1/k 2 ~ 1/uj to the integrand. Each gradient operator 
in ^4(i.2) leads to a factor of k x . Expanding the prop- 
agators to order Q and integrating the resulting term 
over frequency and wavenumber then in general leads to 
an f2( d-2 )/ 2 term in (i-dimensions. Detailed calculations 
confirm this structural argument, see below. The nonan- 
alytic frequency dependence of the FMGM contributions 
to the conductivity is thus as strong as that of the dif- 
fuson contributions. However, we find that the prefactor 
of the FMGM nonanalyticity is of 0(1) in d = 2 + e, in 
contrast to that of the diffuson nonanalyticity, which is of 
0(1/ e). While the two sets of soft modes thus give com- 
parable contributions to the nonanalyticity in the metal- 
lic phase in d — 3, the FMGMs do not lead to a In SI term 
in d = 2. As in the case of the DOS, the conclusion is 
that the FMGMs are irrelevant for the description of the 
metal-insulator transition in d = 2 + e. 

For the technical derivation of these res ults, w e start 
by explicitly writing out ^4(i,2) given by Eq. (4.10a). Note 
that according to the argument given above, the singular 
terms arise from the wavenumber and frequency depen- 
dencies of the soft modes, and all other k and u> depen- 
dence can be neglected. Using this we can localize ^4(i,2) 
in space and time and obtain, 



E 



.4' 



(1,2) m Jm 



(4.13) 



with, 



niri2 rs,ij k cr— ± {3 

x;(<5A)^,„ i+ „ 2 (k)i(<5A)^ + „ 2 , li+ffm (-k) . (4.14a) 
Here k = (k x , k y , fc z ), and N is the matrix element, 

rs N n= E ^ (T a T a T b T r T c T s ) tr (SaSbStScS-j) 
a, b. c— 0,3 

x JdxdyG a n (x)d Xl G b n (x-y)yiG c n (y) , (4.14b) 



with G° and G 3 f rom Eq. fe.3a ), and t± from Eq. ( 4.6e ). 
From Eqs. (4.11) - ( 4.14a| ), and using the factorization 
property of the Gaussian action, we obtain 



iT 



r (i,2) 



mnvm 2 



rs-" m r's' n[ y2 



E E E E-' 

1 ,n 2 ,n' 1 ,n' 2 r,s,r',s' a, a' 

o' 1 



E kx k ' x E 

/3f3> 



k.k 



<;(^)^„ 1+n2 (k);',(5A)«f;, +n ,(kO> G 



x<i(^ + „ 2 ^ i+ _(-k)f,(M)^ +n ,„,_ ff , m (-k')) 
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+<;(^A)^„ 1+n2 (k)f,( ( 5A)f;; + „,^_ (T , m (-k')) G 

x < l(SAti: +nMm (~V l, (SA)% n , +n , (k')> c 

(4.15) 

Note that Eqs. ( t4.14b|) and ( |4.15| ) show that the FMGMs 
and diffusons are not coupled together: If i equals one or 
two, then j must be one or two; similarly, if i is zero or 
three, then j must be zero or three for a nonzero ,J 'JV. 
In what follows, we restrict ourselves to the FMGM con- 
t ribut ions, so we set i,j = (1, 2) and (3 = (3' = a in Eq. 
(4.15). For later reference we also note the symmetry 
relations 



12 jya 

> rs n 



11 N a = 22 N a 

rs n rs n 

33 JV ?i — 00 Jv n> j 30 Jv n — 03 iv ri 



21 jya 

rs n i 



(4.16a) 
(4.16b) 



Next we separate the <5A flucuations into soft SQ fluc- 
tuations and massive <5A fluctuations using Eq. (3.4). To 
make our procedure more transparent, we will initially 
integrate out the massive modes in saddle-point approx- 
imation, i.e. w e simpl y drop the 5A. We will improve on 
this in Section |IVB 3 below. In addition to this approxi- 
mation, we keep only terms that contribute to nonanalyt- 
icities of 0(ft( d - 2 )/ 2 ) in cr( O). T his procedure amounts 
to the replacement, in Eq. ( 4.15 ) and for i, j = 1, 2, 



v N 4 y N I B 

rs ly n rs JV nl Dn ' 



SA => SQ 



(4.17a) 



with 



from Eq. ( 3.8h| ). In terms of Q-fluctations, the leading 
FMGM contribution to <7(i,2) thus is 



CT(1,2) 



iT 



flVm 2 



tj TV 

. rs A ni 



E E E E-' 

,r',s' 



(B ni ) 2 {B^fv 2 



k,k' 



< r mZn l+ n 2 (k) {5Q)Zn' 1+ n> ( k ')) G 

x ( i(SQ)n 1 +n 2 ,n 1 +o-m(~ k) ° s , (SQ) n ^ +n i^ n ^ _ a i m (— k )) G 

+<;w)r,« 1+ « 2 (k) j ;w)^x 1 -.'™(- k ')) G 

X (i(3Q)Z+n 2 ,n 1 +crm(-k) \, (<5Q)"" n , +n , (k')) G . 

(4.18) 

We define 

<;w)r„ +m (k) i{sQ)«?n- +m -V)) G = (27r) d 5(k+k') 

xC^(n,7z + m,k|7i',n' + m') . (4.19a) 



Notice that we have not yet implemented the frequency 
restrictions that we used to define our Gaussian theory 
in Sec. Ill, so the Clj{n-\ , n 2 , k |n,cj, n 4 ) are related to the 
propagators of Eq. ( 3.7a ) by 



C^(™i,™2,k|n3,n 4 ) = — 
lb 



6(ni - n 2 ) Q(n 3 - n 4 ) 

. . — 1 aa.aa 
X J A/I 

rs^ l nin2,n 3 ri4 



(k) 



+ 



x6(ni - n 2 ) 6(n 4 - n 3 ) (-) s 
x8(n 2 - n x ) 6(n 3 - n 4 ) (-) r 
xe(n 2 -n 1 )e(n 4 -n 3 )(-) H 



ij » > — 1 aa,aa , , 
rs'' l nin2 ,n 4 ii 3 I / 



^ A/T 1 aQ ' Qa (V) 
I rs J l n 2 n 1 ,n 3 n i \ I 

/ i 

+ \ / + \ 



jj . j — 1 aa.aa , . 
X rs- M n2ni,n4n 3 i K J 



(4.19b) 



Here we have use d the symmetry relations of the Q- 
matrices, E qs. ( 2.5). We further define the quantity M% bc 
by (cf. Eq. (|4.14b| )) 



M« 6c - | dx dy (x) X1 G* (x - y) yi G c n (y) . 



(4.20) 



Carrying out th e spin a nd qu aterni on su ms in Eq. (4.18) 
and using Eqs. ( |4.16| ), ( fll9| ), and ( |4.20| ), one obtains 



B n = ^- n (k = 0) = i^[^(p)-^(p)] , (4.17b) ^ 



-16T x ^ 
Q,m 2 2— 1 



E 



y ^ s 2 s 2 , 

k ™1 ™i CTCT' 



E™' 



(o-Cqo - ^33) (ct'Cqo - <JC\1 

-{vCoo - 0-^33 ) (cr'C 30 + o-C 03 ) 



33 i 



(4.21a) 

where the frequency labels on the C are C(ni,ni + 
n 2 , k|n' 1; nj + n' 2 ) for the first bracket in each of the four 
terms, and C(n\ + n 2 , n\ + am, — k|"i + «' 2 , «i — cr'm) for 
the second bracket, and 



AT(D 



M 333 _ M 300 

M° 33 + M 303 



M° 30 
M 330 



M° 03 
M„ 000 



(4.21b) 
(4.21c) 



In writing Eqs. (4.21a) we have rearranged the frequency 
labels and made approximations that do not affect the 
leading nonanalytic terms we are aiming to calculate. In 
particular, we have dropped dependences of the prefac- 
tors B n and Nn ' on the "hydrodynamic" frequency la- 
bels n 2 and m. 
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We next perform the sums over a and a' in Eq. (4.21a). 
It turns out that none of the terms that involve G rl or 
C 21 contribute to the leading nonanalytic frequency de- 
pendence of the conductivity. We are thus left with only 
one term, which depends on ^Mr 1 only. Keeping again 
only leading terms, the latter can be expressed in terms 
of the fundamental Goldstone propagator 



g(k,iQ n ) =g + (k,n n )+g-(k,Q n ) 
-2dcp 2 
~ ^ n + (cp) 2 ' 



(4.22) 



with g ± from Eq. fl3.12a| ). Using the Eqs. (§/7|) - ( |3j| ) 
in Eq. (4.21), all of the resulting terms are proportional 
to constants times integrals over products of Goldstone 
propagators. We find 



fT(i ^ )(f7)= i3l /( ^ m ^ o+z0) 

where 

W - ^J2 k " r^e(n 2 ) 5 (k,o; n2 ) 

k n 2 

x g(k,uj n2 +n m ) + g(k,uj n2 



(4.23a) 



and 



a = Tj2£L(k = 0)N^/Bt 



(4.23b) 



(4.23c) 



which is k-independent in our saddle-point approxima- 
tion. Such a k-dependence indeed arises from our for- 
malism if we keep the massive, <5A-fiuctuations that we ne- 
glected so far for simplicityE 



Consider Eq. ( 4.6b ) again. 
By k eepi ng 6 A in the decomposition of A into A and Q, 
Eq. (3.4), and lumping it into the Green function, we 
obtain, upon integrating over 5A in G aussi an approxi- 
mation, a generalization of Q ± in Eq. (4.24) which has 
a k-dependent magnetic self-energy. This procedure still 
constitutes an approximation, as determining the pref- 
actor exactly would require an exact treatment of the 
massive modes. It serves to demonstrate, however, that 
the prefactor is in general nonzero, and our initial null 
result indeed derives from an oversimplifying approxima- 
tion. To avoid misunderstandings, we also point out that 
the leading frequency dependence we determine is exact, 
and only the prefactor we cannot calculate exactly. 

The above procedure produces a wavenumber depen- 
dent self-energy, and hence a nonzero a, even for our 
model with a point-l ike spin-triplet interaction ampli- 
tude r'*', Eq. (2.10c). A wavenumber dependent 
would of course also lead to a wavenumber dependence 
of the self-energy, but capturing this effect would require 
a generalization of our matrix fields. We conclude that 
the prefactor a in Eq. (4.23a) is in general nonzero and 
nonuniversal, as it depends on microscopic details like 
the precise strucure of the interaction amplitude. 



4- The FMGM-induced nonanalyticity 



Notice that S^ n (k = 0) = S nn (k = 0), which follows from 
Eqs. (fT~ 



Simple asymptotic analysis shows that /(£!) — 1(0) oc 
f2 d / 2 . However, before we discuss th e integ ral in detail, 
the prefactor, denoted by a in Eq. ( 4.23a ), warrants a 
closer look. 



3. The prefactor of the nonanalyticity 



We now discuss the prefactor a defined in Eq. (4.23c). 
To leading order in a disorder expansion, it suffices to 
calculate a in the clean limit. Simple manipulations yield 



a 9 51 y a2 



(g+(k)) 2 k x ^-g-(k) 



(&„-(k)) 2 fc,^-e+(k) 



(4.24) 



In the saddle-point approximation for the <5„, Eq. 



( 2.13c ), that we have employed so far, one finds a = 0. 
However, this is an artifact of that approximation. To 
see this, we first point out that a necessary and sufficient 
condition for o ^ is a wavenumber dependence of th e 
'magnetic' piece of the self energy, A„ in Eq. (2.13c), 



We fi nally n eed to consider the integral I(iQ, m ) defined 
in Eq. ( 4.23b ). At T = 0, the relevant dimensionless 
integral is 



J(fi) = 



dk k 



d+l 



u>* - 
k 2 



k 4 



(uj + fl) 2 + k 4 



(w - fl) 2 + k 4 
(4.25a) 



First consider the integral in d = 2. Subtracting the 
value at O = 0, it is easy to see that the leading frequency 
dependence is linear, i.e. there is no term oc f21nf2. More 
generally, standard asymptotic analysis yields 



J(fl -> 0) = J(0) 



2 3+d/2 S i n ^d/A) 



(4.25b) 



Consistent with the absence of a logarithm in d = 2, the 
prefactor of the fi rf / 2 nonanalyticity is finite in the limit 
d — > 2. The frequency sca le for the nonanalyticity is given 
by the quantity A, Eq. ( 2.15| ), which is proportional to 
the magnetization. 

At nonzero temperature, one finds the same qualita- 
tive behavior, since the Goldstone modes, which are the 
source of the nonanalyticity, have the same functional 
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form at all values of T, as long as one stays in the mag- 
netic phase £3 This is in sharp contrast to the diffuson- 
induced weak-localization nonanalyticities, which are cut 
off by a nonzero temperature since the diffusons acquire 
a mass at T > 0. Collecting everything, we obtain our 
final result, 



Rea(Q -» 0) = cr B 



1 



Cd_ 
e F T 



(Q/A) 



(d-2)/2 



(4.26) 



which is valid for 2 < d < 4. C d > in Eq. ([Of]) is 
a positive constant (for fixed dimensionality d) which is 
nonuniversal; it depends, for instance, on the wavenum- 
ber dependence of the spin-triplet interaction amplitude 
as was discussed in the preceding subsection. The d- 
dependence of Cd is nonsingular, in particular, C2 is a 
finite number. 



V. DISCUSSION 

In this paper we have developed a general theoretical 
framework for disordered itinerant ferromagnets. While 
microscopic details like band structure etc. could be 
built into the theory (they would enter in the form of 
more complicated Green functions), our method is par- 
ticularly well suited for studying universal properties that 
are due to the soft modes in the system and that are in- 
dependent of the details on microscopic scales. We have 
therefore restricted ourselves to one parabolic band, and 
have studied the influence of the soft modes on the trans- 
port and thermodynamic properties. Our most impor- 
tant conclusion is that, in disordered systems with ferro- 
magnetic long-range order, there are two distinct families 
of soft modes that contribute to the leading nonanalyt- 
icities that lead to generic scale invariance. One family 
consists of the diffusive "diffusons" that also exist in the 
absence of magnetic long-range order, while the other 
are the magnetic Goldstone modes or spin waves that 
are characteristic of magnetically ordered systems. Even 
though, there are crucial differences between the effects 
of the two types of soft modes, which we now discuss in 
some detail. 



As we have seen in Sec. IV, the one-loop correction 
to the DOS (Fig. |](a)) and the simple loop contribution 
to the one-loop correction to the Boltzmann conductiv- 
ity (the first diagram in Fig. |2|) are simply a frequency- 
momentum integral over a soft propagator. For interact- 
ing electrons, the inter actin g part of the diffuson prop- 
agator, i.e. E in Eq. ( 3.8a ) in the diffuson channels, is 
essentially a diffusion propagator squared, i.e., it scales 
like 1/lu 2 ~ 1/k 4 . The integration then leads to a term 
proportional to S! (d ~ 2) / 2 in 2 < d < 4, and InO in_4 = 2. 
This is the usual weak-localization contributionO The 
FMGM contribution, on the other hand, consists of a 
frequency-momentum integral over a single propagator 
that scales like l/u> ~ 1/k 2 , and is hence less singular. 



As a result, there is no FMGM contribution to the lead- 
ing nonanalyticity in the DOS. 

For the "football" contribution to the conductivity 
(the second diagram in Fig. |J) the situation is differ- 
ent, since both internal propagators can be Goldstone 
modes. Our explicit calculation has shown that this term 
indeed contributes to the leading, 0(fl^ d ~ 2 ^ 2 ), nonana- 
lytic frequency dependence of the conductivity. In d = 3, 
the contributions from the diffusons and the FMGMs are 
thus qualitatively the same. However, in d = 2 the dif- 
fusons lead to the well-known "weak-localization" loga- 
rithmic frequency dependence, while the FMGM contri- 
bution does not lead to a logarithm. A related feature is 
that the prefactor of the nonanalyticity in d = 2 + e goes 
like 1/e for the diffuson-induced term, but is of 0(1) in 
the FMGM case. 

An interesting technical feature is the fact that in the 
case of the diffusons, both the simple loop and the "foot- 
ball" contribute an 0(Sl^ d_2 ^ 2 ) dependence of the con- 
ductivity, while for the FMGM only the latter does. As 
mentioned above, the different results from the simple 
loop can be understood by realizing that the interact- 
ing part of the diffuson is really a diffusion propagator 
squared that scales like 1/k 4 , while the FMGM scales 
like 1/k 2 . A question that arises then is why the dif- 
fusons in the "football" do not yield a much stronger 
singularity, of 0(fl( d ~ 6 ^ 2 ). The technical reason is that 
the diffusons come with a dditio nal frequency restrictions, 
as was discussed in Sec. IIIB. These frequency restric- 



tions lead to additional frequency factors in the inte- 
gral that determines the "football" contribution, while 
no such singularity-protecting effect occurs in either the 
simple loop for the diffusons or in any of the FMGM 
contributions. 

Another very important difference between diffusons 
and FMGMs is that the former are soft only at T = 0, 
while the latter remain soft at all T below the Curie 
temperature. As a result, the nonanalytic frequency de- 
pendence coming from the diffusons gets replaced by an 
analogous nonanalytic temperature dependence at T > 0. 
The one induced by the FMGM, .an the other hand, is 
still of 0{^ d - 2 ^ 2 ) even at T > OB 

The sign and strength of the nonanalytic term in Eq. 
( 4.26[ ) is interesting and can be understood on very gen- 
eral grounds. First we remember that this nonanalyticity 
in frequency also exists at finite temperature. This im- 
plies that its origin is classical in nature. The sign of 
classical "mode coupling" or "generic scale invariance" 
contributions to the transport coefficients is determined 
by the structure of the hydrodynamic equations describ-, 
ing the long wavelength dynamics of the soft modes.cfJ 
For the spin density dynamics relevant here, the leading 
nonlinear coupling in the long wavelength limit contains 
two gradients, as it does-in classical dissipative systems 
such as Lorentz models.Ej As in the Lorentz models, the 
soft modes therefore have a localizing effect, i.e. the con- 
ductivity decreases with decreasing frequency. In con- 
trast, in classical fluids the leading nonlinear coupling 
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contains only one gradient, and the-prefactor of the non- 
analyticity has a delocalizi ng sig n.EII We further note that 
the nonanalyticity in Eq. ( 4.26| ) is stronger than the one 
in classical Lorentz gases by a factor of ft. It is, in fact, 
as strong as the generic scale invariance effects in clas- 
sical fluids, even though the structures of the respective 
nonlinear couplings suggest just the opposite. The reso- 
lution of this apparent paradox lies in the fact that there 
is long-range order in our system, which is manifested by 
the Goldstone modes. In a mode-coupling calculation, 
the transverse spin suspectibilty, \t 1/k 2 , appears as 
a multiplicative factor in the integrands and effectively 
leads to the same power of the frequency as in the fluid 
case. It is interesting to note that a similar mechanism 
in certain liquid crystals leads to mode coupling effects 
that are even stronger than in classical fluids .E3 For these 
systems the standard hydrodynamic description breaks 
down for all d < 5. Here, the analogous critical dimen- 
sion is d = 2, see Eq. ( 4.26| ). 

There are some important experimental consequences 
of these results. One is that, as far as the weak- 
localization properties in d — 2 are concerned, one does 
not expect any differences between metals with long- 
range ferromagnetic order and nonmagnetic metals in 
an external magnetic field. The reason is that the only 
difference in the soft mode structures of the two sys- 
tems are the Goldstone modes, but they do not con- 
tribute to the leading nonanalyticities in d = 2. Indeed, 
transport measurements on thin ferromagnetic films have 
found just the ordinary weak-localization effects due to 
the diffusons, with_no additional effects from the Gold- 
stone modes at all.Ej Our theory explains this null result, 
which a priori is very surprising. For bulk materials, on 
the other hand, our theory predicts that the FMGMs do 
contribute a term that is qualitatively the same as the 
weak-localization effect, and furthermore that this term 
will remain a nonanalytic frequency dependence even at 
T > 0. 

Another consequence is the suggestion that the 
FMGMs will not be important for the critical behavior 
at the metal-insulator transition from a FM metal to a 
FM insulator, at least close to two-dimensions. The rea- 
son for this is the absence of a diverging prefactor of the 
frequency nonanalyticity as d — > 2. This divergence is 
known to drive the transition in d = 2 + e, and since the 
FMGMs do not contribute to it one expects the univer- 
sality class of the transition to be unchanged by the pres- 
ence of ferromagnetic long-range order SB We will inves- 
tigate this point in (II) and find that the metal-insulator 
transition on the background of ferromagnetism is indeed 
closely related to the one of nonmagnetic electrons in an 
external magnetic field. 
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